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ABSTRACT 

For the quadruple gravitational lens PG 1115+080, we combine recent 
measurements of the time delays with new lens models to determine the Hubble 
constant Hq. We explore the effects of systematic uncertainties in the lens 
models on the estimates of Hq, and we discuss how the uncertainties can be 
reduced by future observations. We find that the lens cannot be fit by an 
isolated lens galaxy, but that it can be well fit by including a perturbation 
from the nearby group of galaxies. To understand the full range of systematic 
uncertainties it is crucial to use an ellipsoidal galaxy and to let the group 
position vary. In this case, the existing constraints cannot break degeneracies 
in the models with respect to the profiles of the galaxy and group and to the 
position of the group. Combining the known time delays with a range of lens 
models incorporating most of the plausible systematic effects yields Hq = 51^11 
km s~^ Mpc~^. The constraints on the lens models, and hence on Hq, can be 
improved by reducing the standard errors in the lens galaxy position from 50 
mas to ~ 10 mas, reducing the uncertainties in the time delays to ~ 0.5 days, 
and constraining the lens mass distribution using HST photometry and the 
fundamental plane. In particular, the time delay ratio r^BC = ^tac / ^^ba may 
provide the best constraint on the mass profile of the galaxy. 



Subject headings: gravitational lensing - cosmology: distance scale - galaxies: 
structure - quasars: individual (PG 1115+080) 
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1. Introduction 

A general consensus is emerging from local studies that the Hubble constant Hq lies in 
the range 65-75 km s~^ Mpc~^ (e.g. Freedman, Madore & Kennicutt 1996), based primarily 
on Cepheid distances to nearby galaxies combined with type la supernovae to reach regions 
of pure Hubble flow (e.g. Riess, Press & Kirshner 1996). The resulting age estimates are 
in weak conflict with the estimated ages of the oldest stars for almost all cosmological 
models (e.g. Bolte & Hogan 1995). Measurements of time delays in gravitational lenses 
can directly determine Ho over cosmological distances (Refsdal 1964), thereby avoiding the 
complicated calibration problems that plague local distance estimates. In a gravitational 
lens, the ray trajectories of the multiple images have different geometric lengths and pass 
through different parts of the gravitational potential, so the light travel time is different for 
each image. The light travel time is inversely proportional to Hq, so combining observed 
time delays between images with a model of the gravitational potential gives the Hubble 
constant. Unfortunately, application of this technique has been slowed by the difficulties of 
determining the time delays and of flnding good models for the gravitational potential. 

The first gravitational lens discovered, 0957+561 (Walsh et al. 1979), was also the 
first lens for which a time delay was measured (e.g. Vanderreist et al. 1992; Lehar et al. 
1992; Press, Rybicki & Hewitt 1992a,b; Schild & Thomson 1995; Pelt et al. 1996; Kundic 
et al. 1996); the most recent measurement yielded the time delay At = 417 ± 3 days at 
95% confidence (Kundic et al. 1996). The system has been modeled extensively (e.g. Young 
et al. 1980; Borgeest & Refsdal 1984; Narasimha, Subramanian & Chitre 1984; Greenfield, 
Roberts & Burke 1985; Falco, Gorenstein & Shapiro 1991; Kochanek 1991b; Bernstein, 
Tyson & Kochanek 1993; Grogin & Narayan 1996). Most recently, Grogin and Narayan 
(1996) considered a spherical softened power-law model for the primary lens galaxy and an 
external shear for the surrounding cluster. The best-fit galaxy model had a dark matter 
halo whose mass increases slightly faster than isothermal, M(r) cx r° with 1.07 < a < 1.18 
at 95% confidence. This model gave Ho = (1 - «;)(At/1.1 yr)"^ km s"^ Mpc-\ 

where k > is an inherent degeneracy due to the mean surface mass density of the cluster 
(Falco, Gorenstein & Shapiro 1985). Kundic et al. (1996) estimated k from the observation 
of the cluster by Fischer et al. (1997) and used the model of Grogin & Narayan to infer 
Ho = Q4:± 13 km s~^ Mpc^^ (95% confidence); independently, Falco et al. (1997) used the 
velocity dispersion of the galaxy to infer iJo = 62 ± 7 km s^^ Mpc^^. We note that Grogin 
& Narayan (1996) briefiy considered an approximate elliptical model of the galaxy, but they 
fixed the ellipticity and position angle to match the observed isophotes and did not explore 
the additional freedom in the models (and hence in Ho) due to treating the galaxy and 
cluster as independent sources of shear. 
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Recently the four-image gravitational lens PG 1115+080 (Weymann et al. 1980; see 
Figure |1|) became the second lens for which time delays were measured. Schechter et al. 
(1997) measured the time delay between images C and B to be Ar^c = 23.7 ±3.4 days and 
the time delay between C and mean of the close pair A = Ai + A2 to be Ar^c = 9.4 ± 3.4 
days, giving a time delay ratio tabc = ^tac / ^tba = 0.7 ± 0.3. Bar-Kana (1997), however, 
reanalyzed the data to show that including the correlations in the photometric errors gives 
Atbc = 25.0l};f days and rABC = l-13to;i7 (statistical), with a ~ 0.2 systematic uncertainty 
in Tabc associated with different assumptions about how to treat the photometric errors. 

PG 1115+080 is a promising candidate for combining time delays with lens models to 
determine Hq, because four-image lens geometries can constrain some aspects of lens models 
better than two-image lens geometries (Kovner 1987; Kochanek 1991a). PG 1115+080 has 
been modeled extensively. The first models demonstrated that simple non-axisymmetric 
galaxy models could qualitatively reproduce the image configuration (Young et al. 1981; 
Narasimha, Subramanian & Chitre 1982). Further models showed that an isothermal sphere 
or a point mass with a variety of quadrupole structures could reproduce the image positions 
to within ground-based observational errors, but that the observations were not accurate 
enough to distinguish between the various monopole and quadrupole forms (Kochanek 
1991a). More accurate observations with the Hubble Space Telescope (Kristian et al. 1993, 
hereafter K93) have provided better observational constraints. Keeton, Kochanek & Seljak 
(1997a) showed that the quasar image positions and fluxes and the galaxy position cannot 
be fit by models with a single shear axis (such as an ellipsoidal galaxy or a circular galaxy 
with an external shear), but can be fit by models allowing two independent shear axes 
(such as an ellipsoidal galaxy with an external shear). Schechter et al. (1997) demonstrated 
that the nearby group of galaxies seen by Young et al. (1981) can provide the required 
external shear. In addition, Schechter et al. (1997) used their measurements of the time 
delay with a simple lens model treating the galaxy and group as singular isothermal spheres 
to infer Hq = 42 km s~^ Mpc~^, with a 14% uncertainty from the time delay and with 
unknown uncertainties from the lens model. Their value for Hq is in strong conflict with 
local estimates of Hq = 65-75 km s^^ Mpc~^, so it is crucial to understand the systematic 
uncertainties in the lens models and their effects on the inferred value of Hq. In addition, 
if we hope to use lensing to measure Hq with precision comparable to the distance ladder, 
then we must understand which future observations can reduce the uncertainties in the lens 
models. 

In this paper we examine the systematic uncertainties in the value of Hq inferred 
from PG 1115+080 by exploring the uncertainties in the lens models. Since present 
observations of the galaxy and group do not directly constrain their mass distributions 
and only weakly constrain their positions, we postulate a wide range of plausible galaxy 
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and group properties. Specifically, we model the galaxy as an ellipsoidal mass distribution 
with a range of profiles including constant mass-to-light ratio models, standard dark matter 
models, and more centrally-concentrated models. We model the group as a circular mass 
distribution with an extended (isothermal) or concentrated (point mass) profile. We then 
study the ability of the data to constrain the models and show that present constraints 
leave significant degeneracies in the models. We examine the implications of the models 
and degeneracies for Hq, and discuss future observations that can break the degeneracies 
and improve the constraints on Hq. In §2 we present the data and methods used in the 
models. In §3 we consider a broad range of models using only the primary lens galaxy and 
show that they cannot fit the data. In §4 we study models adding the nearby group and 
show that a variety of galaxy+group models give good fits. We use a Bayesian analysis 
with the most physically plausible models to quantify the present systematic uncertainties 
in Hq. In §5 we illustrate an independent way to break one degeneracy by using stellar 
dynamics and the fundamental plane of elliptical galaxies to identify which galaxy profiles 
are physically plausible. In §6 we summarize our results for Hq as well as the prospects 
for breaking the degeneracies through better observations of the quasar fluxes, the galaxy 
position and photometry, and the time delays. 

2. Data and Methods 

PG 1115+080 consists of four point images of a 2;^ = 1.722 quasar with separations 
~ 2" surrounding a galaxy at redshift 0.295 ± 0.005 (Weymann et al. 1980; K93; Angonin- 
Willaime et al. 1993). There is a nearby group of galaxies (Young et al. 1981) at redshift 
0.304 (Henry & Heasley 1986). A schematic diagram of the lens and the galaxies is shown 
in Figure |l[ For reference, in an f2o = 1 cosmology 1" at the group redshift is 2.77h^^ kpc. 
Other cosmological scale factors used in the lensing analysis are given in Table |1] for several 
cosmologies. Throughout the text we use an Qq = 1 cosmology. 

Computing Ho from a gravitational lens requires time delays and a lens model. The 
time delays were taken from Bar-Kana's (1997) analysis of the data of Schechter et al. 
(1997). The B — C time delay is the longest and best-resolved time delay. Bar-Kana (1997) 
found Atbc- = 25.0ti;7 days, consistent with the value 23.7 it 3.4 days from Schechter et al. 
(1997), and this result depended only weakly on assumptions about the photometric errors. 
By contrast, the other time delays and hence the time delay ratio tabc = ^tac I ^'^ba are 
less well determined. Bar-Kana (1997) found tabc = l-13loj7, conflict with the value 
0.7 ± 0.3 from Schechter et al. (1997), and this result did depend on assumptions about 
the photometric errors. A range of assumptions indicated that Bar-Kana's (1997) result 
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for tabc "was more robust than the result of Schechter et al. (1997), and that systematic 
uncertainties in r^BC were at the level of ~ 0.2 (Bar-Kana 1997). 

To determine the lens model, we fitted lens mass distributions to reproduce the quasar 
images positions and fluxes and the galaxy position from the Hubble Space Telescope 
observations by K93 (see Table §). The relative coordinates have an uncertainty of 5 mas 
for the quasar images and 50 mas for the lens galaxy. The quasar flux ratios are less well 
determined because of the source variability, microlensing, and extinction. To account for 
this, we broadened the flux error bars to 20% to encompass the range of observed variability 
(see Keeton & Kochanek (1996) for a summary). We also considered the effects on our 
conclusions of making the flux error bars smaller (see §6). K93 estimated the upper limit 
on the flux of a faint central image to be 1-2% of the flux of the brightest image, and we 
included this constraint by setting the limit on the flux of the central image to be ± 2% of 
the flux of the brightest image. 

Time delays offer an independent constraint on lens models, but only if more than 
one is known. The first time delay is used to determine Hq, and the rest are combined 
into iJo-independent ratios that constrain the models. However, because of the systematic 
uncertainties in the time delays for PG 1115+080 we used them only to determine Hq, 
not to constrain the models. For each model, we used the B — C time delay to compute 
i^O; which we express as Hq = 100/iBC'(ArBc/25.0 days)~^ km s~^ Mpc~^. We used the 
other known time delay in the ratio tabc only as a qualitative check of the consistency of 
the models. We did, though, examine the effects of including the tabc constraint in the 
Bayesian analysis of Hq in §4.2. 

To determine whether a lens mass distribution is consistent with the data, we studied 
the images it can produce. The lensing properties of a surface mass distribution S are 
described by its lensing potential ip determined from the two-dimensional Poisson equation 
V^V = 2S/Ecr (e.g. Schneider, Ehlers & Falco 1992), where the critical surface density for 
lensing in angular units is 
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rn = c/Hq is the Hubble radius, and Di, Ds, and Dis are angular diameter distances to 
the lens, to the source, and between the lens and the source, respectively. Values of the 
distance ratio are given in Table |l|. The potential ip describes the mapping between the 
source and image planes through the lens equation 

u = X - Vip{x), (2) 

where x is an angular position in the image plane and u is an angular position in the 
source plane. A source at u maps to images at Xj that are roots of the lens equation. The 
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images are deformed by the magnification tensor Mij where Mj~^ = Sij — d'^ip/dxidxj. 
Because surface brightness is conserved, the total magnification of an image is the ratio of 
the area of the magnified image to the area of the source; so the magnification factor is 
M = det(Mjj). The ray trajectories of the images have different geometric lengths and pass 
through different parts of the gravitational potential, so the light travel time is different for 
each image; the time at image position x is 



r(f) 



1 + zi DiD, 



D, 



X — u 



— %Ij{x) 



(3) 



and the time delay between an image at Xi and an image at Xj is Ar^j = r(xj 
Note that r factors into a piece that depends on the lens model times a piece that depends 
on cosmology and scales as Hq^ (see Table |l| for values in different cosmologies). If we let 
D = D/th he 3ja. angular diameter distance scaled by th (and hence independent of i^o), 
then we can use eq. (|]) to find Hq from a lens model and time delay. 
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(4) 



From Table |l|, changing the cosmology can change the inferred value of Hq by up to 7%. 
We modeled the primary lens galaxy with ellipsoidal surface mass densities of the form 



where 



m 



l + ecos2(e-^. 



(5) 



e is a natural ellipticity parameter such that the axis ratio is (1 — e)^/^/(l + e)^^"^, and 9^ 
is the orientation angle of the major axis. We will quote 6^ as a standard position angle 
measured North through East. We considered two classes of models. First, we used the de 
Vaucouleurs (1948) model as the prototypical constant mass-to-light ratio (M/L) model for 
early-type galaxies, with surface mass density 



b exp |-/c (m^/r^)^/^| 
dvv exp{— fcf^/^} 



(6) 



where k = 7.67, is the effective (or half-light) radius, and b is the deflection scale such 
that the total mass is 
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for b and r^. in arcseconds. Second, we used softened power-law models with surface mass 
density 



[s^ + m 



2U-a/2' 
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where a is the power law exponent such that M(r) oc asymptotically, s is a core radius, 
and h is the deflection scale. For a = 1 it is a softened isothermal model, whose lensing 
properties have been studied by Kassiola and Kovner (1993) and Kormann, Schneider & 
Bartelmann (1994). For a = it is a modified Hubble model (Binney & Tremaine 1987). 
We also studied the more centrally-concentrated a = — 1 model. For a = 1, 0, —1 we found 
analytic expressions for the deflection and the magnification, and for other values of a we 
used numerical integrals. 

We modeled the group as a single halo in which the observed galaxies are embedded. 
For simplicity we used only circular mass distributions. We considered a standard dark 
matter model treating the group as a singular isothermal sphere (SIS), and we studied the 
effects of making it more centrally concentrated by considering the limit of a point mass. 
The lensing potential for a point mass at the origin is ?/^(r) = }? Inr. 

For a given lens mass distribution, we solved the lens equation (^ to map a source at 
u to its images at Xj. We varied the position and flux of the source and the parameters 
of the mass distribution using the "amoeba" downhill simplex method (Press et al. 1992) 
to minimize the residuals in the image plane. We included in the the constraints from 
the quasar positions and fluxes and the galaxy position. The galaxy position and the four 
quasar positions and fluxes provide 14 constraints. A galaxy model with re or (s, a) fixed 
has five free parameters [xgai-, e, 6*^, and 6). A group model has one parameter {bgrpj if 
the position is fixed and three parameters if the position is variable. The source has three 
parameters {u and the flux). Thus a lens model using an isolated lens galaxy with fixed 
or (s, a) has N^oj = 6 and a lens model using the galaxy plus a movable group has Ndof = 3. 

3. Results: Isolated Galaxy 

We first considered simple lens models using only the primary lens galaxy. Previous 
studies showed that the K93 data cannot be fitted by an isolated lens galaxy if the galaxy is 
represented by a point mass or a singular isothermal mass distribution (Keeton et al. 1997; 
Schechter et al. 1997). To see whether some other isolated galaxy could fit the data, we 
considered a range of softened power-law models, as well as de Vaucouleurs constant M/ L 
models. For all models we used an ellipsoidal mass distribution as the source of shear. 

For the de Vaucouleurs models, there are no observational estimates of the effective 
radius Vg, so we tabulated results as a function of Ve (see Figure 0). To estimate the 
plausible range of r^, we note that an galaxy has an effective radius of (4 ± l)h^^ kpc 
(Kormendy & Djorgovski 1989; Rix 1991), which corresponds to 1'.'4 at the redshift of 
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PG 1115+080. In addition, in §5 we use stellar dynamical arguments to estimate that 
lies in the range I'.'O <^re <^ 4'.'0. With this range of effective radii, no de Vaucouleurs model 
gives an acceptable fit. The x^/^dof is 68 at = I'.'S (4.2/i~^ kpc) and is even larger at 
smaller effective radii. Although x^/^dof decreases for larger r^, it is still 18 at Vg = 10" 
{28h~^ kpc). All physical parameters inferred from the models are unreasonable: for 
Te = 0'.'5-4'.'0, the models require a very flattened galaxy {b/a ~ 0.1-0.4) and imply a huge 
Hubble constant {hsc ~ 2.3-1.3). Moreover, the ratio of the model Ai — C and B — Ai 
time delays is tabc ~ 5.4-2.0, compared with 0.7 ± 0.3 for the time delays of Schechter et 
al. (1997) or 1.13 ± 0.2 for the time delays of Bar-Kana (1997). A simple constant M/L 
galaxy cannot describe PG 1115+080. 

For the softened power-law models, we tabulated the results in the plane of the core 
radius s and the power-law exponent a (see Figure 0). Generally the softened power- law 
profiles fit much better than the de Vaucouleurs profile but still do not give a good fit. The 
best-fit model has l^dof = 11-7, is nearly isothermal (a = 1.002), and has a large core 
radius {s = 0'.'45 = 1.2h"^ kpc). This model implies an axis ratio b/a = 0.75, a Hubble 
constant hsc = 0.50, and a time delay ratio vabc = 1-44 that is marginally consistent 
with Bar-Kana (1997) given his systematic uncertainties. As s or a decreases, the models 
become more centrally-concentrated and imply larger values for Hq; this is because Hq 
scales as 1 — where is the convergence at the critical radius (Falco et al. 1985), and the 
models with steeper profiles have a smaller kq. The models are not tightly constrained, with 
the la region allowing core radii 0'.'28 (0.8/i~^kpc) <^ s <^ 0'.'55 (1.5/i^^ kpc), and power-law 
exponents 0.7 ^ « ^ 1.6. Profiles steeper than the modified Hubble profile {a = 0) are 
ruled out at better than 99.9% confidence. Over the la region, tabc is well constrained 
and varies by only ~ 0.1, but hsc is poorly constrained and varies from 0.2 to 0.7. For 
reference, the singular isothermal ellipsoid model has x^ l^dof = 25, h/a = 0.50, hsc = 0.94, 
and Tabc = 1-85. 

It is surprising that the models require such a large core radius even for the isothermal 
profile. We expect a small core radius for two reasons. First, observations of galactic 
cores show that the luminosity densities of elliptical galaxies do not have fiat cores but 
instead have central cusps (Gebhardt et al. 1996). Second, it is generally argued that small 
core radii are needed to fit lens data, both because almost all known lenses lack a central 
image, and because stellar dynamics, lens statistics, and other lens models are consistent 
with isothermal mass distributions with small core radii (Narasimha, Subramanian & 
Chitre 1986; Kassiola & Kovner 1993; Wallington & Narayan 1993; Kochanek 1993, 1995, 
1996; Grogin & Narayan 1996). Contrary to these expectations, PG 1115+080 cannot 
be fit by an isolated lens galaxy with a small core radius. It is difficult, though, to draw 
strong conclusions from this result because the "best-fit" model is still not a good fit 
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{x^/Ndof — 11-7), while once we include the perturbation from the group we will find good 
fits (x^/Noof ^ 1) with galaxies that do have a small core radius. 

We conclude that PG 1115+080 cannot be fit by an isolated lens galaxy, and thus it is 
not reasonable to use such a model to estimate Hq. 

4. Results: GalEixy + Group 

Keeton et al. (1997a) found that while PG 1115+080 could not be fit by an isolated 
isothermal galaxy, it could be well fit by an ellipsoidal galaxy with an independent external 
shear. Schechter et al. (1997) pointed out that the group seen by Young et al. (1981) was 
correctly located to be the source of the shear, and by modeling the galaxy and the group 
as two singular isothermal spheres they inferred Hq — 42 km s~^ Mpc~^. This value of 
Hq is unexpectedly low, but it is difficult to judge its significance because Schechter et al. 
(1997) did not explore the systematic uncertainties in the lens models and their effects 
on the inferred Hq. In order to understand the systematic effects, we now examine the 
galaxy+group models in detail. We explore the freedom in the models (and hence in Hq) 
and discuss the prospects for better constraining the models with future observations. First 
we consider models with the group position fixed to examine the effects of the profiles of the 
galaxy and the group, and then we consider models in which the group position is allowed 
to vary. 

4.1. Fixed group 

The group mass distribution is described by a position and any parameters associated 
with the radial profile. Let dgrp be the distance from the primary lens galaxy to the 
group, and bgrp be the critical radius of the group. To lowest order the group is a 
perturbation that can be characterized by its convergence Kgrp and shear 7^,^^ together 
with weaker nonlinear terms (Falco et al. 1985). For a singular isothermal sphere (SIS) 
group -fgrp = Kgrp = bgrp/2dgrp, while for a point mass group 7^^^ = b^^p/dlrp i^grp = 0- 
Since ^grp scales with bgrp/ dgrp, we expect that requiring an external shear 'jgrp ~ 10% 
(Kochanek 1991a; Keeton et al. 1997a; Schechter et al. 1997) will produce (to lowest order) 
a degeneracy between the group's position and its mass. The non-linear terms may break 
the degeneracy (see §4.2), but for simplicity we first consider models with the group position 
fixed. In this way we can focus on the effects of different profiles for the galaxy and the 
group. 
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The three group galaxies and the lens galaxy are at essentially the same redshift and 
are probably physically related, so we placed the group at the flux-weighted centroid of all 
four galaxies (C4 in Figure 0). The Gunn r magnitudes of the three group galaxies were 
taken from Young et al. (1981); the r magnitude of the primary lens galaxy was estimated 
from the F785LP magnitude (K93) by using the galaxy evolution models of Bruzual Sz 
Chariot (1993) to compute the color of an E/SO galaxy at a redshift of 0.3 (see Keeton, 
Kochanek & Falco 1997b for details). We found the flux centroid to be dgrp = 14'.'5 from 
the primary lens galaxy at a position angle 6grp = — 117?1 (North through East). 

Figure | shows the results for de Vaucouleurs, softened isothermal (a = 1), modified 
Hubble (a = 0), and a = —1 models of the primary lens galaxy with either an SIS or a 
point mass group, and Table ^ summarizes the ^ind the physical parameters for the best 
fits. All eight classes of models give good fits, with x^/^dof < 1 for all but the model with 
an a = — 1 galaxy and a point mass group. Since no model with an isolated galaxy could 
do better than x^/^dof ~ 10, while numerous models with the galaxy supplemented by the 
group produce /Ndof ^ 1, we confirm the results of Keeton et al. (1997a) and Schechter 
et al. (1997) that including the group is crucial to obtaining a good fit. We note that the 
best model of Schechter et al. (1997), which treated the galaxy and group and singular 
isothermal spheres, gave x^ l^dof = 5. Our models give x^ l^dof < 1 despite the addition 
of the flux constraints because we allowed the galaxy position to vary, and because we 
included a second shear by allowing the galaxy to have an ellipticity. As we noted in Keeton 
et al. (1997a), the key to fitting many of the quadruple lenses is having two independent 
shear axes. 

The values for Hbc implied by these models vary significantly with the galaxy and 
group profiles (see Table Qualitatively, the variation makes sense: Hq scales with 1 — 
(Falco et al. 1985), and the convergence kq at the critical radius decreases as the galaxy 
or group becomes centrally concentrated. Unfortunately, the data cannot constrain the 
profiles for two reasons. First, although the four-image configuration constrains the total 
mass within the critical radius, it does not constrain the distribution of the mass (Kochanek 
1991a; Wambsganss & Paczyhski 1994). Second, although the quasar images constrain 
the shear from the group, they do not constrain its mass or convergence. As a result, we 
can have a wide range of galaxy and group profiles that are all consistent at the la level, 
and by using different profiles we can produce hsc anywhere from 0.4 to more than 0.8 
without significantly changing the goodness of fit. As we show in §5, though, not all of the 
profiles corresponds to galaxies that are physically plausible. The centrally concentrated 
models that allow the higher values for Hq correspond to galaxies that do not lie on the 
fundamental plane of elliptical galaxies. Restricting attention to plausible dark matter and 
constant M/L models restricts the range of hsc to 0.4 <^ Hbc ^ 0.7. 
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4.2. Movable group 

While fixing the group at the flux centroid C4 allowed us to isolate the effects of the 
profiles of the galaxy and group, there is no strong observational evidence requiring it. 
Hence the group position is an additional systematic uncertainty whose effects on Hq must 
be examined. We studied the effects of the group position for two galaxy models: a dark 
matter (isothermal) model and a constant M/L model. For the constant M/L model we 
used a modified Hubble profile because, although it is not as good a representation of galaxy 
luminosity profiles as the de Vaucouleurs profile, it does have an analytic deflection formula. 
We neglected the more centrally concentrated galaxy models because, as we show in §5, 
they require galaxies that are unphysical. For the group we again used the SIS and point 
mass models. Physically we expect the group to be described by an isothermal dark matter 
halo, but by including the point mass group we can examine the effects of making the 
group more concentrated. Since these two galaxy profiles and two group profiles span the 
expected range of mass profiles, they should span the expected range of Hq. We examined 
all four combinations of profiles; representative results are shown in Figures § and |^, where 
Figure |^ shows the extended profiles (isothermal galaxy and SIS group), and Figure || shows 
the concentrated profiles (modified Hubble galaxy and point mass group). 

The singular isothermal ellipsoid (dark matter) galaxy permits good fits for a wide 
range of group positions (see Figure With an SIS group, the best-fit model is at 
dgrp = 25'.'2 and 9grp = —125°, has = l-^T for Ndof = 3, and implies hBc = 0.47. 
But at the Icr level, Og^p can range from —150° to —110° and dg^p can be as small as a 
few arcseconds. The Hubble constant varies considerably over the allowed region, from 
hsc = 0.2 for distant groups with 9grp ~ —113°, to hsc = 0.7 for nearby groups or for 
distant groups with 9grp ^ —145°. 

Note that there is no formal upper limit on dgrp from the x^; because for a distant 
group the non-linear terms are weak and the perturbation is equivalent to a simple 
external shear, which we know fits well (Keeton et al. 1997a). However, the mass of the 
group increases with its distance {Mgrp oc bgrpdgrp ~ d'^j,p), so requiring that the group 
have a reasonable mass can constrain its position. It is convenient to study the mass of 
the group in terms of its mass-to-light ratio (M/L), and it is convenient to express the 
M/L of the group in units of the M/L of the galaxy because this ratio is independent of 
assumptions about the cosmology and about the K and evolutionary corrections to the 
luminosities. The models with the group in a band from the centroid C4 to the best-fit 
group position have 5 ^ {M/L)grp/{M/L)gai ^ 25. The galaxy mass within the critical 
radius (radius tq = 1'.'14 = 3.16h~^ kpc) is 1.22 x lO^^h'^ Mq, with a 5% uncertainty 
due to the position of the group, a 10% uncertainty (upwards) due to making the group 
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more centrally concentrated, and a 7% uncertainty (upwards) due to the cosmological 
model. Using the K93 F785LP central aperture magnitude together with color, K, and 
evolutionary corrections from the galaxy evolution models of Bruzual & Chariot (1993), we 
estimated the present B magnitude of the galaxy to be —19.2 mag for an f2o = 1 cosmology 
(see Keeton et al. 1997b for details), although the effective aperture of the K93 magnitude 
is unclear. This gives {M/L)gai — 18, which is consistent with other lens models (Keeton 
et al. 1997b) and lens statistics (Kochanek 1993, 1996; Maoz & Rix 1993). This in turn 
implies 90 ^ {M/L)grp ^ 450, which is consistent with observed group mass-to-light ratios, 
e.g. 150 ^ M/L ^ 350 (Ramella, Pisani & Geller 1997). It is difficult at present to apply 
strong constraints from this reasoning, but it is clear that the group positions that give 
good fits also give reasonable group masses and that the group should not be much closer 
or much farther away. 

Figure ^ shows that at present the constraints on Hq are relatively weak. Nevertheless, 
the models imply several physical properties that may be better constrained by future 
observations. First, the galaxy positions differ at the level of tens of mas. While the K93 
position is well within the range of good fits, the 50 mas error bars are too broad to rule out 
many models. Reducing the error bars on the galaxy position to 10 mas or better would 
greatly improve the constraints on the models. Second, the galaxy axis ratios vary. The 
galaxy is strictly circular if the group is at dgrp = 13'.'5, 9grp = —115° (which reproduces 
the double SIS model of Schechter et al. 1997). For other group positions the galaxy's 
ellipticity and orientation adjust so the combined shear from the galaxy and group gives a 
good fit. The ellipticity needs to vary by only a few times 0.05 in order to produce good 
fits for a wide variety of group positions. Finally, the time delays vary. Although it was not 
included as a constraint, the ratio tabc is marginally consistent with the value 1.13 ± 0.2 of 
Bar-Kana (1997) given his systematic uncertainties of ~ 0.2. Note that no models produce 
a time delay ratio consistent with the value tabc = 0.7 it 0.3 of Schechter et al. (1997). The 
time delay ratio could provide a strong constraint on the models if the systematic effects in 
the observed ratio were understood and if the uncertainties were reduced to ~ 0.05. This 
would require roughly a factor of four reduction in the uncertainties in the time delays, or 
uncertainties of ^ 0.5 days. We note also that the variation of the A2 — Ai time delay with 
group position differs from that of the A — C and B — A time delays, so this time delay 
could provide an independent constraint. However, Ar^^^^ = 5 ± 1 hours over the la region 
so this time delay would need to be known at the level of ~ 0.5 hour in order to constrain 
the models. 

The results are similar for the modified Hubble model (constant M/L) galaxy, except 
that good fits are limited to a narrower range of group positions (see Figure With a 
fixed core radius s = 0'.'2 {0.55h~^ kpc) and a point mass group, the best-fit model is at 
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dgrp = 28'.'3 and 9grp = —118°, has = 1-70 for Ndof = 3, and implies hsc = 0.67. At 
the la level, 9grp is restricted to —120° ^ 9grp ^ —112° and the group can be no closer 
than dgrp = 11". As with the isothermal galaxy, the Hubble constant varies considerably 
over the la region, ranging from 0.3 to 0.9. However, again physical parameters including 
the galaxy position, the galaxy axis ratio, and the time delays may be constrained by 
further observations to improve the constraints on Hq. In particular, note that the modified 
Hubble galaxy predicts a time delay ratio that is closer to the observed value of Bar-Kana 
(1997) than the ratio predicted by the isothermal galaxy. Although it is not shown, the 
predicted time delay ratio is essentially independent of the profile of the group. Thus if 
Bar-Kana's (1997) value remains valid as its uncertainties are reduced, then it may provide 
the strongest probe of the galaxy profile. We will return to this point below. 

Given the wide range of group positions that provide good fits, it is difficult to use 
the statistic to place limits on Hq. We can, however, use a Bayesian analysis to give a 
reasonable estimate of the systematic uncertainties. Using Bayes's theorem we can convert 
the probability of the data given the parameters (given by e~^^^^) into the probability of 
the parameters given the data. We can then compute the probability distribution for Hq by 
integrating over the group position, weighted by a "prior" probability distribution which we 
took to be a circular Gaussian distribution centered on the flux centroid C4 with standard 
deviation given by the rms distance of the four galaxies from the centroid. This is roughly 
the same as the range of group positions permitted by the group mass-to-light ratio. With 
the probability distribution for Hq we can characterize the systematic uncertainties in Hq 
due to the group position, and also estimate the relative likelihoods of the four classes of 
galaxy/group profiles. In addition, because the four classes of profiles span the range of 
physically plausible models (a dark matter or constant M/L galaxy, and an extended or 
concentrated group), we can combine them to produce a "total" probability distribution 
that includes the systematic effects of both the group position and the galaxy/group 
profiles. Figure |^ shows the normalized probability distributions and inferred values of Hq 
for the four classes of models, together with the total probability distribution. Since the 
time delay ratio r^BC provide a strong constraint on the galaxy profile but is still 
relatively uncertain, we have computed the probability distributions with and without the 
formal vabc constraint from Bar-Kana (1997). 

The Bayesian analysis emphasizes two important features of the uncertainties in Ho 
due to systematic uncertainties in the lens models. First, using two independent shears (an 
ellipsoidal galaxy and a movable group) strongly affects Hq. Dark matter models that have 
only a single variable shear, such as the models of Schechter et al. (1997) with a circular 
galaxy or the models of §4.1 with the group position fixed, give Ho ~ 40 km s~^ Mpc~^. 
By contrast, dark matter models with two shears give Hq = 58t^\l km s~^ Mpc~^. Thus to 
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understand the systematic uncertainties in Hq it is important to consider models with two 
independent shears. Moreover, since lens models that require two independent shears often 
have a degeneracy between the shears (Keeton et al. 1997a), such models will generically 
produce large uncertainties in Hq that can be reduced only by improving the constraints on 
the models. 

Second, the time delay ratio Vabc is an important constraint on the profile of the 
galaxy. The current estimate reduces the range of models that are consistent with the data 
(especially the models that predict large Ho), and reduces the probability of the isothermal 
galaxy compared with the Hubble galaxy. Without the tabc constraint, the isothermal 
galaxy is more likely than the Hubble galaxy by a ratio of 4 : 1, largely because the 
isothermal galaxy allows such a wide range of group positions. By contrast, with the tabc 
constraint the Hubble galaxy increases in likelihood because it produces values for r^BC 
that are more consistent with Bar-Kana's (1997) value. Reducing the uncertainties in r^BC 
thus should help discriminate between the dark matter and constant M/L galaxy models 
and provide the best probe of the galaxy profile. Using the formal Tabc constraint from 
Bar-Kana's (1997) present estimate, the total probability distribution gives Hq = km 
s~^ Mpc~^, where these error bars incorporate most of the systematic uncertainties in the 
lens models due to the group position and the galaxy/group profiles. 



5. Stellar Dynamics 

To this point we have evaluated models purely on their ability to fit the gravitational 
lens data, and we have found that the data leave significant degeneracies in the models. One 
degeneracy is related to the galaxy profile, but here we can apply independent constraints 
to try to break the degeneracy. Specifically, we can consider whether the galaxies required 
to fit the lensing data are consistent with stellar dynamics. For example, we suspect that 
the centrally concentrated models — the models that allow high values for Hq — may be 
unphysical. 

The absence of information on the optical structure of the lens galaxy limits the 
conclusions that can be drawn from stellar dynamical models because such models require 
the luminosity distribution to estimate the stellar velocities. We assume that the galaxy 
is an early-type galaxy because of its high mass and the expected dominance of lens 
statistics by early-type galaxies. Early-type galaxies are known to obey a strong correlation 
between velocity dispersion, effective radius, and magnitude known as the fundamental 
plane (Djorgovski & Davis 1987; Dressier et al. 1987), and we can estimate whether the 
lens galaxy as constrained by the lens models can lie on the fundamental plane. Kochanek 
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(1993, 1996), Breimer & Sanders (1993), and Grogin & Narayan (1996) have previously 
explored using stellar dynamical models as an added check on lens models. 

For our comparison sample we adopted the data of J0rgensen, Franx & Kjaergaard 
(1995a,b, hereafter JFK). Figure |^ shows the JFK sample in the space of logcr and log re 
shifted to the redshift of PG 1115+080. The velocity dispersion a was estimated for a fixed 
metric aperture of 3'.'4 in diameter at Coma (corresponding to 0'.'4 at PG 1115+080 for 
fio = !)• In order to compare the K93 and JFK photometry, we converted the JFK Gunn 
r magnitudes to the K93 F785LP band by using the galaxy evolution models of Bruzual 
& Chariot (1993) to compute the K and evolutionary corrections and the r — F785LP 
color (see Keeton et al. 1997b for details). We assumed that the luminosity density can be 
modeled by a Hernquist (1990) distribution with Hernquist scale length a = 0.45re, and 
we calculated only isotropic stellar dynamical models. From the models in §4.2, we know 
that with an SIS group the mass inside the ring of images (radius tq = r.'14 = 3.16h^^ 
kpc) is 1.22 X lO^^h^^ Mq, with a 5% uncertainty due to the position of the group, a 10% 
uncertainty (upwards) due to making the group more centrally concentrated, and a 7% 
uncertainty (upwards) due to the cosmological model. 

Figure [^superimposes the predicted PG 1115+080 aperture velocity dispersion as a 
function of effective radius on the JFK sample. The self-gravitating Hernquist (constant 
M / L) models have diverging central velocity dispersions for small effective radii because the 
mass is fixed by the lens model, forcing a oc for effective radii smaller than the ring 
defined by the images. The effective radius is restricted to the rough range I'.'O ^ re ^ 4'.'0 
if the galaxy is to lie on the fundamental plane; then from §4.1 the de Vaucouleurs models 
with fixed group give 0.5 ^ hsc ^ 0.7 with an SIS group and 0.6 ^ hsc ^ 0.9 with a 
point mass group. Surprisingly, the dark matter models tend to have dispersion estimates 
that move along the fundamental plane. Nonetheless, many of the centrally-concentrated 
lens models that lead to high values for Hq require a galaxy off the fundamental plane; for 
example, the a = —1 model with s = O'.'l that produces hBc = 0.77 lies to the right of the 
JFK sample. These fundamental plane constraints can be made quantitative once we know 
the lens galaxy structure and effective radius. 

Note that the Hernquist models require a total lens galaxy magnitude of /(F785LP) ~ 17 
in order to lie on the fundamental plane, whereas K93 estimated a central aperture magnitude 
for the lens galaxy of J(F785LP) = 18.36 mag. We experimented with the K93 images 
and found that, due to the wings of the original WFPC point-spread function and to the 
huge contrast between the quasar and galaxy surface brightnesses, the data are consistent 
with lens models having /(F785LP) ~ 17 and — l'.'5-2'.'0. In addition, Keeton et al. 
(1997b) have remarked that the K93 magnitude estimate gives a luminosity well below 
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that expected from the "Faber-Jackson" type relation between image separation and lens 
luminosity that other lenses obey. A higher luminosity for the PG 1115+080 lens galaxy 
would move it closer to the trend followed by other lenses. Improved HST photometry will 
reveal whether this is indeed the case. 

Although this analysis is limited at present by the lack of optical data on the lens 
galaxy, it does reveal that the centrally-concentrated galaxy models that give high values for 
Hq are physically implausible. The dark matter models are consistent with stellar dynamics, 
as are the constant M/ L models if the effective radius is in the range 1"0 ^ Te ^ 4'.'0. 



6. Discussion 

By combining a lens model of PG 1115+080 with the observed time delays of Bar-Kana 
(1997), we can infer a value for Hq that is independent of the standard distance ladder. The 
resulting value depends on the lens model, so we explored a range of models to understand 
the systematic uncertainties in the models, their effects on Hq, and the types of future 
observations that can reduce the uncertainties. We found that PG 1115+080 cannot be fit 
by an isolated lens galaxy, but that it can be well fit ix^ /N^of < 1) by including the effects 
of the nearby group. To understand the full range of model uncertainties, it is important 
to use an ellipsoidal galaxy and a movable group. Since the present observational data 
do not constrain the galaxy and group profiles or the group position, we studied various 
profiles and group positions and found good fits with Hq ranging from roughly 30 to 90 
km s~^ Mpc~^. However, many of these models are unacceptable on physical grounds. For 
example, some of the models that give high values for Hq require a galaxy mass profile that 
is more centrally concentrated than typical luminosity profiles and correspond to a galaxy 
that does not lie on the fundamental plane. If we restrict the models to a plausible constant 
M/L model and a dark matter model, we can use a Bayesian analysis to characterize the 
systematic uncertainties in Hq relating to degeneracies in the models. We find Hq = 
km s~^ Mpc~^, where the error bars incorporate the uncertainty in the measured time delays 
as well as the uncertainties in the lens models due to the group position and the galaxy 
and group profiles. They do not, however, take into account the systematic uncertainties 
in the value of the time delay ratio tabc from Bar-Kana (1997). In addition to the formal 
uncertainties, there may be a 5-10% uncertainty due to mass fluctuations from large-scale 
structure (Bar-Kana 1996; Wambsganss et al. 1997), and a 7% uncertainty (upwards) due 
to the cosmological model (Table |I]). 

The uncertainties in Hq can be reduced with better observational constraints, but 
perhaps surprisingly not with better relative quasar positions and fluxes. In most models 
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the quasar positions are already overfit, and an order of magnitude improvement in the 
uncertainties (to 0.5 mas) is not practical. The quasar fluxes, by contrast, are dominated by 
systematic effects from microlensing and extinction rather than by measurement errors. We 
accounted for these effects by using flux error bars of 20%; merely reducing the flux error 
bars by a factor of two changes the absolute but has little effect on the bcst-flt model 
parameters. Improving the constraints from the quasar fluxes would require understanding 
the systematic effects of microlensing and extinction, and in particular understanding why 
the K93 value of the A2/A1 flux ratio (0.66 in V and 0.70 in I) differs from the best-fit 
models (typically 0.90) and from theoretical expectations of a value near unity. One way 
to avoid problems with microlensing and extinction would be to measure radio fiux ratios; 
unfortunately PG 1115+080 is not a strong radio source (flux < 1.5 mJy; Weymann et al. 
1980). 

The most promise for improving the constraints comes from the galaxy position and 
the time delays (see §4.2). Reducing the error bars on the galaxy position from 50 mas to 
10 mas or better will rule out many models, and such a reduction should be possible with 
new WFPC2 observations. Reducing the error bars on the time delays by a factor of four 
(to ~ 0.5 day in the time delays or to ~ 0.05 in tabc) will provide a useful independent 
constraint. In fact, the time delay ratio vabc may provide the best probe of the galaxy 
proflle. The effects of the galaxy position and the time delay ratio can be seen by noting 
that, together with the convergence at the ring of images, they account for almost all of 
the variation in Hq from model to model. If we combine the models from §4.1 (different 
galaxy proflles with varying core radius and flxed group) and the models from §4.2 (singular 
isothermal ellipsoid or modifled Hubble model galaxy, with movable group), then the 
implied Hubble constant is correlated with the galaxy position and the time delay ratio, 

Hbc = 0.83 + O.OOBAxgai - O.OOeAygai + 0.10(rABC - 1-4) (1 - ko) (9) 

where Axgai and Aygai are the galaxy position in mas relative to that of K93, and 
'^o = i^gai + i^grp IS the total convergcncc at the ring of images due to the galaxy and the 
group. For a singular isothermal galaxy Kgai = 1/2 at the ring of the images; for an SIS 
group Kgrp = jgrp ~ 10% in PG 1115+080, while for a point mass group Kgrp = 0. This 
correlation is meant only as a qualitative guide to the quantities that will best constrain 
the Hubble constant, and it does not address the degree to which the convergence Kq can 
be determined from improved data. 

If even with improved constraints we still cannot determine the galaxy and group 
profiles and the group position directly, then perhaps we can impose external considerations. 
First, as in §5 we can use stellar dynamics and the fundamental plane to identify 
which models are physically plausible. At present this analysis can only rule out very 
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concentrated models, but with luminosity profile data for the galaxy it may provide a more 
quantitative constraint. Unfortunately, directly measuring the central velocity dispersion 
in PG 1115+080 will be very difficult due to the high quasar/galaxy contrast. Second, we 
could require consistency with other gravitational lens statistics and models (Grogin & 
Narayan 1996; Kassiola & Kovner 1993; Kochanek 1993, 1995, 1996, WaUington & Narayan 
1993; Maoz & Rix 1993) and with observations of elliptical galaxies (Fabbiano 1989; Rix 
1996) to say that lens galaxies have significant dark matter and small core radii. In this 
case, the models from §4.2 give i^o = 44 ± 11 km s^^ Mpc~^ using the tabc constraint. 
Finally, there are now two lenses (0957+561 and PG 1115+080) for which time delays have 
been measured. Both systems exhibit a degeneracy between the mass distribution and the 
Hubble constant (e.g. Grogin & Narayan 1996), but by requiring the lenses to agree on both 
the Hubble constant and the typical mass distributions of galaxies we may be able to break 
the degeneracies. 
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Table 1. Cosmological scale factors 
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(0.4,0.6) 


(0.2,0.8) 


1" 

DiDs 

(l+2i) DiD, 
c As 

y 
y 


{h-^ kpc) 

days arcsec"^) 
{W^h-^MQ arcsec-2) 
(/i g cm~^) 


2.771 
0.139 
30.49 
0.326 
0.888 


2.968 
0.149 
32.73 
0.350 
0.831 


3.080 
0.145 
31.68 
0.339 
0.747 


3.219 
0.145 
31.80 
0.340 
0.686 



Note. — Quantities are computed for Zg = 1.722, zi = 0.304. The first two cosmological 
models are FRW cosmologies with the specified fig. The last two cosmological models are 
flat cosmologies with + Aq = 1- 



Table 2. Observational data 





x{") 


y (") 


F555W (mag) 


F785LP (mag) 


Al 


-1.294 


-2.036 


16.90 


16.12 


A2 


-1.448 


-1.582 


17.35 


16.51 


B 


0.362 


-1.949 


18.87 


18.08 


C 


0.000 


0.000 


18.37 


17.58 


G 


-0.355 


-1.322 




18.36 



Note. — Data from HST WFPC observations by Kristian ct al. (1993). x is approximately 
west, y is approximately north. The internal position uncertainties arc 5 mas for the quasar 
images and 50 mas for the galaxy. Formally, the relative fluxes are uncertain by 1.5% 
in I (F785LP) and 3% in V (F555W), while the zero-point for magnitudes is uncertain 
by 0.3 mag. The positions and Gunn r magnitudes of the group galaxies are as follows 
(Young et al. 1981): Gl= {dgrp = 23'.'5,9grp = -118°, r = 18.96), G2= (12'.'0, -95°, 20.04), 
G3= (18'.'9, -131°, 20.53). 
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Table 3. Galaxy + fixed group models 
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Note. — Results from the best-fit galaxy + fixed group models. The group position is 
fixed at the flux- weighted centroid of the four galaxies (C4 in Figure 1). The de Vaucouleurs 
models were computed only for discrete spaced by 0'.'5. The softened power-law models 
were computed with s varying continuously. The Xpos x'jiux indicates the contribution 
to the from the quasar positions and fluxes, respectively. Each model with fixed Vg or s 
has Nfiof — 5. We do not give error bars because they would not encompass the effects of 
allowing the group to move. 
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Fig. 1. — Schematic diagram of the gravitational lens PG 1115+080 and the nearby group. 
The crosses represent the four quasar images (K93). The circles represent the galaxies (Young 
et al. 1981), with the areas of the circles indicating the relative fluxes. The plusses represent 
( "C3" ) the flux- weighted centroid of the three group galaxies and ( "C4" ) the flux- weighted 
centroid of all four galaxies. The physical scale indicated is for an = 1 cosmology. 



-25- 




Fig. 2. — Results for models using only the primary lens galaxy. The galaxy is represented 
by an ellipsoidal mass distribution with the specified density profile. Each model with fixed 
Te or (s, a) has N^of = 6. Left: de Vaucouleurs profile. The heavy solid line is the 
using the axis scale on the left. The light solid fine is the Hubble constant expressed as 
Hq — 100/iBc(^T-Bc/25.0 days)~^ km s~^ Mpc~^, and the dashed fine is the axis ratio of the 
galaxy; both use the axis scale on the right. Right: Softened power-law profile. The solid 
lines are contours of drawn at Ax^ = 2.30, 4.61, 6.17, 9.21, 11.8, and 18.4, the Icr, 90%, 
2(7, 99%, 3(J, and 99.99% confidence levels for two parameters. The dotted lines are contours 
of hsc- The best-fit model (marked with a triangle) is at s = 0'.'449, a = 1.002 and has 
X^ = 46.9 and hec = 0.50. 
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Fig. 3. — Results for models using the primary lens galaxy plus the nearby group. The 
group position is fixed at the fiux- weighted average of the four galaxies (C4 in Figure 1). 
The primary lens galaxy is represented by an ellipsoidal mass distribution and the group 
is a represented by a circular mass distribution. Top: the group is represented by an SIS. 
Bottom: the group is represented by a point mass. The heavy solid line is the using the 
axis scale on the left; each model with fixed Te or s has N^of = 5. The filled and open points 
indicate the best fit and the A^^ = 1 and A^^ = 4 limits. The light solid line is hsc ^-nd 
the dashed line is the axis ratio of the galaxy, both using the axis scale on the right. 



-27- 



Galaxy = Sinigular Isothermal Ellipsoid, Group = SIS 
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Fig. 4. — Results for models representing the primary lens galaxy as a singular isothermal 
ellipsoid and the group as an SIS. With the group position fixed a model has N^of = 5. (a) 
The sohd lines are contours of drawn at Ax^ = 0.2, 0.4, . . . (light) and at Ax^ = 2.30, 
4.61, 6.17, 9.21, 11.8, and 18.4 (heavy), the Icr, 90%, 2(t, 99%, 3(7, and 99.99% confidence 
levels for two parameters. The dashed fines are contours of Hbc- The positions of the group 
galaxies and of the fiux-weighted centroid (C4 in Figure 1) are indicated, (b) Contours of the 
time delay ratio tabc- The heavy contours indicate the la range l.lSlo^y from Bar-Kana 
(1997). (c) Contours of the distance A^qj = \xmod — Xobs\ (in mas). The heavy contour 
indicates the range allowed by the Kristian et al. (1993) error bars, /^yai < a/2 x 50 mas. (d) 
Contours of the galaxy ellipticity 1 — b/a. In the upper right corner the galaxy is oriented 
approximately North-South, and everywhere else it is approximately East-West. The best- 
fit model (marked with a triangle) is at dgrp — 25'.'2 and 9grp = —125° and has x^ — 1-77, 
hBC — 0.47, Tabc — 1-50, Agai — 25 mas, and b/a — 0.85. 
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Galaxy = Modified Hubble (s = 0.2), Group = Point Mass 
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Fig. 5. — Results for models representing the primary lens galaxy as a modified Hubble 
model with core radius s = 0'.'2 (0.55/i~^ kpc) and the group as a point mass. The panels are 
the same as in Figure 4, except that the range of position angles is smaller. Again in (d) the 
galaxy is oriented approximately East-West in the lower left and North-South in the upper 
right. The best-fit model (marked with a triangle) is at dgrp = 28'.'3 and 6grp = —118° and 
has = 1-70, Hbc = 0.67, tabc = 1-43, Agai = 16 mas, and b/a = 0.92. 
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Fig. 6. — Normalized probability distributions for Hq. The distributions were computed 
from a Bayesian analysis of the four classes of models discussed in §4.2. In the top panel the 
Bayesian analysis does not use the constraint from tabc, and in the bottom panel it does 
use tabc- The relative probabilities of the four models and their implied values for Hq are 
given in the key. The spread within each distribution is due to the 6% uncertainty in the 
observed time delay and to the degeneracy in the group position. The spread between the 
distributions is due to the degeneracies in the profiles of the galaxy and the group. 
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Fig. 7. — Stellar dynamical limits. The points show the effective radii and central dispersions 
of the JFK sample of cluster galaxies, with the effective radius rescaled to the redshift of 
PG 1115+080. The heavy solid lines show the central dispersion for the lens galaxy assuming 
the same metric aperture as a function of rg; for the softened power-law models, the core 
radius s is given in arcsec. The light solid lines show the F785LP magnitude estimates 
for the lens galaxy from the fundamental plane. The K93 value /(F785LP) = 18.4 is a 
central aperture magnitude and represents a lower bound. We found that the K93 images 
are consistent with a galaxy having rg ~ I'.'S and /(F785LP) ~ 17. 



